Polaronic behavior of undoped high-T c cuprates 
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We present angle-resolved photoemission spectroscopy (ARPES) data on undoped La2Cu04, indi- 
cating polaronic coupling between bosons and charge carriers. Using a shell model, we calculate the 
electron-phonon coupling and find that it is strong enough to give polarons. We develop an efficient 
method for calculating ARPES spectra in undoped systems. Using the calculated couplings, we find 
the width of the phonon side band in good agreement with experiment. We analyze reasons for the 
observed dependence of the width on the binding energy. 



An important issue in the discussion of high-T c 
cuprates is the relative importance of couplings to 
phonons and spin fluctuations. Recent angle-resolved 
photoemission spectroscopy (ARPES) work for undoped 
or weakly doped Ca2-zNa :r Cu02Cl2 shows polaronic be- 
havior y|. A broad k-dependent boson side band was 
observed, while a quasi-particle, if present, was too weak 
to be seen. This raises important questions about the 
nature of the bosons coupling to the electrons. Here 
we present experimental data for La2- x Sr a; Cu04, showing 
similar boson side bands. Using a shell model, we calcu- 
late the electron-phonon coupling for undoped La2Cu04. 
We find that the coupling is sufficiently strong to give 
polaronic behavior. We develop a method for calculat- 
ing spectra of a t-J model with many phonon modes and 
at finite temperature. The calculated couplings give side 
band widths in good agreement with experiment. The de- 
pendence of the width on the binding energy is analyzed. 
This work shows that the electron-phonon coupling plays 
an important role for properties of undoped cuprates. 

The photoemission measurements were carried out on 
beamline 10.0.1 at the Advanced Light Source. The pho- 
ton energy was 55 eV and the energy resolution was 
w 18 meV. The La2Cu04 single crystals were grown by 
the traveling solvent floating zone method. The samples 
were cleaved in situ in vacuum with a base pressure bet- 
ter than 4x 10 -11 Torr and measured at a temperature of 
ss 20 K. The data are shown in Fig. ^ There are broad 
boson side bands at about 0.5 eV binding energy, similar 
to the results for Ca 2 - x Na a; Cu02Cl2 0]. 

To see if these boson side bands can be explained 
by the coupling to phonons, we calculate the electron- 
phonon interaction. We use a shell model Q that suc- 
cessfully describes phonon dispersions in several cuprates 
to obtain the phonon modes in La2Cu04. Within this 
model the change in the electrostatic potential induced 
by each mode is calculated to first order in the atomic 
displacements. This gives rise to a linear coupling be- 
tween doped charges and the phonons. We use the t-J 
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FIG. 1: (a) Photoemission spectra of La2Cu04 along the 
(0,0)-(7r,7r) nodal direction in the first Brillouin zone. The cor- 
responding momentum runs from (0.32,0.32)-7r to (0.76,0.76)7r, 
as indicated by the arrow, (b) To highlight the momentum 
dependence, a "background" given by a spectrum near (ir,n) 
was subtracted from the spectra. The bold curves correspond 
to a momentum near (ir/2,ir/2). 



model to describe the electronic degrees of freedom in 
the CuC>2 planes. Doped holes form Zhang-Rice singlets 
with the spins of Cu d holes with their additional charge 
distributed symmetrically on the four neighboring O p 
orbitals. This charge distribution couples to the modula- 
tion of the electrostatic potential caused by the phonons. 
Since the undoped system is an insulator, screening in 
the shell model is only provided by the displacement of 
the shells against the atomic cores. In addition, there 
is coupling from the modulation of the p-d hopping t p d 
and the charge-transfer energy in the three-band model 
from which the t-J model is derived Q. We include the 
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FIG. 2: Spectral distribution of coupling to singlets (full line) 
and to holes in non-bonding O p orbitals (dashed line). A 
Gaussian broadening (FWHM=3 meV) was used. 

dominant on-site terms. The modulation of the charge- 
transfer energy is calculated from the modulation of p 
and d level energies within the shell model instead of as- 
suming a nearest-neighbor Coulomb repulsion U p d as in 
Ref. |3|. The shell model calculations are done on a (30) 3 
mesh in q space using the high-temperature tetragonal 
structure from Ref. |2| but similar results are obtained 
for the low-temperature orthorhombic structure. The re- 
sults depend on the eigenvectors of the shell model. Since 
the shell model @ describes neutron scattering intensi- 
ties well, the eigenvectors are believed to be accurate . 
The electron-phonon interaction is then given by 

H ep = -L^ 5q „(l - n^V^vQqve^ 11 *. (1) 

This is an on-site coupling to empty sites representing 
singlets in the t-J model. It is linear in the generalized 
phonon coordinates Q ql , and the strength is given by the 
coupling constants g^ v - A phonon mode with eigenfre- 
quency o> q „ is labeled by its wave vector q and branch 
index v. n.$ measures the occupancy at site R^. The 
number of sites is N. 

Introducing the dimensionless coupling constant A = 
2 E q Jffq^l7(8^q^), we obtain A = 1.2 | for t = 
0.4 eV. The criterion for polaronic behavior in the 
Holstein-t-J model is A > A c = 0.4 5]. If the next- 
nearest neighbor hopping integral t' is taken into account, 
A c should increase. The coupling derived for undoped 
La2CuC>4 should, nevertheless, be strong enough to cause 
polaronic behavior, as observed experimentally. 

Fig. |21 shows the spectral distribution 7(w) = 
J2q V \9^y\ 2 / {^uN)5{uj - w q „) of our coupling (full line). 
A peak around 20 meV can be attributed mainly to 
modes involving La whereas the spectral weight at GO- 
TO meV is mainly due to vibrations of the apical O. 
A peak around 85 meV is caused by the planar O 
(half-)breathing mode coupling predominantly via the 



t p d modulation. This energy dependence of the coupling 
agrees qualitatively with the observation of fine struc- 
tures in the electron self-energy of La 2 - 2; Sr ;r Cu04 0, ex- 
cept that the calculation does not give appreciable cou- 
pling around 40 meV. This could be due to surface effects 
or distortions due to doping 0, |j| . Such effects would 
be missing in our calculation, which is performed for an 
ideal La2CuC>4 structure. For comparison, the dashed 
line in Fig. [21 shows the spectral distribution of the cou- 
pling to holes in individual non-bonding O p orbitals. 

To discuss the q dependence of the coupling we sum 
|<?qv| 2 over all modes and the phonon momentum q z per- 
pendicular to the CuC>2 planes. The coupling increases 
with decreasing | (q x , q y )\ although some individual modes 
show different behavior, e.g. the coupling to the O 
(half-)breathing mode peaks around (n, n). 

Based on Ref. we have developed an efficient 
method for calculating ARPES spectra in undoped sys- 
tems at a finite temperature T. Using an adiabatic ap- 
proximation, the spectral function for creating a hole 
with momentum k in an undoped system is written as 

A k (w, T) — JdQ \MQ, r)| 2 Pk(Q, «). (2) 

The integration is over all phonon coordinates Q ql / sum- 
marized in the vector Q. pk(Q,w) is the spectrum cal- 
culated for the electronic system in a distorted lattice 
(characterized by Q). No dynamic phonons are included 
in the calculation of pu(Q, w), since the distortions Q are 
treated as c-numbers and only give rise to static on-site 
energies in the Hamiltonian 0J. 

w , (aT) i._n^«p( 5 3£ I ) (3 ) 

is the squared ground-state wavefunction for 
non-interacting phonons at T = (or the 
corresponding ensemble average for finite T, 
n>qv{T) — (exp(a; qz , / '(k sT)) — is the Bose-Einstein 
distribution). Our approximation neglects the kinetic 
energy of the phonons in the final states and assumes 
that kgT is small compared with the electronic energy 
scale. We use the fact that there is no electron-phonon 
coupling for the initial state in the undoped t-J model. 

In practice, Eq. can be evaluated efficiently by 
Monte Carlo sampling over the phonon coordinates us- 
ing \(f>o(Q,T)\ 2 as the weighting function. Each sampled 
Pk(Q,w) requires only the calculation of a purely elec- 
tronic problem. There is no blow-up of the Hilbert space 
through the inclusion of dynamic phonons that usually 
prevents calculations for strong coupling. 

To obtain ARPES spectra for La2Cu04 we describe 
the electronic system by the two-dimensional t-J model 
choosing t — 0.4 eV and J = 0.3t. pk(Q,u) is calcu- 
lated using exact diagonalization on a 4x4 cluster with 
periodic boundary conditions. To check our method we 
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FIG. 3: ARPES spectra for the undoped system at T = 
for different k normalized to the height of the phonon side 
band. The lower abscissa shows binding energies (BE) and the 
upper abscissa the energies of the final states corresponding 
to the spectral features. The inset shows the dependence of 
the width of the phonon side band on its binding energy. The 
width of the (0, 0) spectrum is poorly defined and not shown. 



have first assumed a simple Holstcin-type of coupling to 
one dispersionless mode as in Ref. |5j and found good 
agreements with their results. But our method allows us 
to include all 21 modes with the q-dependent coupling 
given by Eq. l]| . 

We have calculated Ay(uj : T = 0) for different values 
of k as shown in the main part of Fig. |3J For each spec- 
trum 50000 samples were used. The poles in the sampled 
spectra were broadened by Gaussians with a FWHM of 
33 meV. The spectra show a broad phonon side band 
which disperses like the quasi-particle in the t-J model 
without phonons, as was found in Ref. [jj. As shown in 
Ref. 01 this can be understood easily within our adi- 
abatic approximation. The true quasi-particle peak is 
almost completely suppressed in weight and dispersion 
and determines the energy zero on the lower abscissa in 
Fig. El The half- width (HWHM) of the side band is deter- 
mined on the low binding energy side. We obtain a width 
of 2xHWHM=0.5 eV for k = (7r/2,7r/2), in good agree- 
ment with the corresponding experimental result 0.48 eV 
(cf. bold curves in Fig.^). 

Width and position of the phonon side band as well 
as the position of the quasi-particle peak are shown in 
Fig. 0] as a function of the relative coupling strength A 
(substituting g qi , — > Ag q; , in Eq. Q) for k = (7r/2,7r/2) 
at T = 0. For small A, the position of the quasi- 
particle peak [l3| starts out at the energy obtained 
in the t-J model without phonons and approaches for 
stronger couplings asymptotically the curve given by 
— ^ 2 |5qi>| 2 / LJ qi^ + const. A fully localized hole ob- 
tains this energy gain from interaction with phonons (the 
additional constant depends on the definition of the en- 
ergy offset). The position of the phonon side band shows 
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FIG. 4: Width and position of the phonon side band and the 
position of the quasi-particle peak as a function of the relative 
coupling strength A for k = (tt/2,tt/2). 



only a weak linear dependence on A. The energy differ- 
ence to the position of the true quasi-particle peak de- 
termines the binding energy. The value of almost 1.2 eV 
at full coupling is larger than the experimental value of 
about 0.5 eV (cf. Fig. QJ. Since the binding energy de- 
pends strongly on doping, a possible explanation for this 
discrepancy could be that the experimental samples were 
actually lightly doped. Our shell model calculation could 
also overestimate the coupling strength. As seen in Fig.0J 
the binding energy of the side band has a much stronger 
dependence on A than the width. For A = 0.8, e.g., 
one would find a width of about 0.4 eV and a binding 
energy of about 0.6 eV, both quantities being in reason- 
able agreement with experimental values. As A oc A 2 , 
this would reduce the dimensionless coupling constant to 
A = 0.75. Such a reduction could be due to a slight 
underestimate of the screening in the shell model or an 
overestimate of the coupling to breathing phonons. 

The inset in Fig.|31shows that the width of the phonon 
side band increases roughly linearly with its binding 
energy. This trend has also been found experimen- 
tally for weakly doped Ca2- K Na 2 ,Cu02Cl2 Q. To un- 
derstand this, we introduce rescaled phonon coordinates 
Rqv = \Av/(2ra qi ,(T) + ijQqv- The weighting func- 
tion in Eq. © then becomes invariant under rotations 
in R space, and Eq. J2J can be expressed as a sam- 
pling over directions for fixed R = |R| followed by an 
integration over R. If the R space has dimension d, 
the i?-integral contains a factor R d ~ 1 e~ R , which peaks 
strongly at R = y/(d— l)/2. Fig. El shows how the con- 
tributions to the spectrum evolve with increasing R up 
to R. Poles in the sampled spectra were broadened by 
Gaussians with a FWHM of 11 meV. In the figure, the 
coupling to the q=(0,0) modes is not included, since it 
contributes the same broadening for all k. We focus on 
peaks to the right, at the lowest binding energies. For 
R = 0, which corresponds to zero electron-phonon cou- 
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FIG. 5: Contributions to the spectra with k = (tt/2, it, 2) 
and k = (0, 7r,2) at T — for varying R. Higher spectral 
weight is represented by darker shading. The abscissa shows 
binding energies (BE). Coupling to q = (0,0) modes is not 
included. 

pling, the spectrum for k = (ir/2, tt/2) has a peak at 
smaller binding energy than for k = (0,7r/2). As R in- 
creases and the electron-phonon coupling is switched on, 
both peaks are broadened. This is due to the repulsion 
from higher states, which increases with R, but is differ- 
ent for different directions of R. For the peak with the 
lowest binding energy this broadening is largest, since all 
the other states repel the corresponding state in the same 
direction. Fig. [3| shows, however, the opposite trend for 
the width of the phonon side band in the final spectra, 
due to the following opposing effect which dominates. For 
R ^ spectral weight also appears at energies with no 
peaks in the R = spectra. For R ^ 0, the system is dis- 
torted and the electronic momentum k is not conserved. 
For instance, in case of k = (0, 7r/2), a peak of increasing 
width and weight appears at the energy of the main peak 
in the k = (tt/2, tt/2) spectrum. At R = R, these side 
bands have more or less merged with the main peak and 
effectively add to its width. For k = (tt/2, tt/2) there 
is no such extra contribution for lower binding energies 
whereas for k = (0,7r/2) the width is increased by one 
side band. There are more and more side bands on the 
low binding energy side of the main peak as the bind- 
ing energy of the main peak increases, which increases 
the width of the resulting broad peak on the low binding 
energy side. 

Finally, we have studied the temperature dependence 
of the width of the phonon side band. It is relatively 
weak. For k = (tt/2, tt/2) it increases from 500 meV 
at T = to 565 meV at T = 200 K and 750 meV at 
T = 400 K. The main contribution to the T dependence 
is due to the modes with an energy of about 20 meV. 

To summarize, we have shown that the polaronic fea- 



tures in ARPES spectra of La2Cu04 and other undoped 
cuprates can be attributed to strong coupling between 
the photo hole and phonons. This electron-phonon cou- 
pling is due to the modulation of the electrostatic poten- 
tial and the singlet energy. We have introduced an effi- 
cient method for obtaining ARPES spectra in undoped 
systems and shown that with the derived coupling there 
is good agreement with experimental results. The de- 
pendence of the width of the broad phonon side band on 
binding energy and temperature has also been discussed. 
Our results show the importance of electron-phonon cou- 
pling for the physics of undoped cuprates. 

We thank T. P. Devereaux and P. Prelovsek for useful 
discussions. 
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